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Abstract 

The collision dynamics of ^^02('^S^) 02(^S^) in the presence of a magnetic field is studied 
within the close-coupling formalism in the range between 10 nK and 50 mK. A recent global ah 
initio potential energy surface (PES) is employed and its effect on the dynamics is analyzed and 
compared with previous calculations where an experimentally derived PES was used [New J. Phys 
11, 055021 (2009)]. In contrast to the results using the older PES, magnetic field dependence 
of the low-field-seeking state in the ultracold regime is characterized by quite a large background 
scattering length, o^g, and, in addition, cross sections exhibit broad and pronounced Feshbach 
resonances. The marked resonance structure is somewhat surprising considering the influence of 
inelastic scattering, but it can be explained by resorting to the analytical van der Waals theory, 
where the short range amplitude of the entrance channel wave function is enhanced by the large 
a\,g. This strong sensitivity to the short range of the ah initio PES persists up to relatively high 
energies (10 mK). After this study and despite quantitative predictions are very difficult, it can be 
concluded that the ratio between elastic and spin relaxation scattering is generally small, except for 
magnetic fields which are either low or close to an asymmetric Fano-type resonance. Some general 
trends found here, such as a large density of quasibound states and a propensity towards large 
scattering lengths, could be also characteristic of other anisotropic molecule-molecule systems. 
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I. INTRODUCTION 



Ultracold molecules play a main role in modern physics due to a large number of promising 
applications in quantum information precision spectroscopy [2] and ultracold chemistryjs- 
5|. Optical lattices of ultracold molecules are predicted to be ideally suited for quantum sim- 



ulation of complex quantum systems 
information storage and processing|9 



and the engineering of new schemes for quantum 
10|. On the other hand, creation of a Bose-Einstein 
condensate (BEG) of molecules may enable studies of Bose-enhanced chemistry llj . In 
the context of these studies, molecules must be confined within a trap. For paramagnetic 
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molecules, a magnetic trap is used since molecules in a low-field-seeking state 
trappable provided that their translational energy is lower than the trap depth 
situation could be achieved by direct cooling methods such as Zeeman slowing [14], optical 
Stark deceleration [l^, single-photon cooling 16 1 or sympathetic coohng[17|]. It might also 
be possible to cool the molecules towards the ultracold regime by evaporative coo! 
As it is well known, this was the successful method for achievement of BEG of atoms 

Molecular collisions are fundamental in this context, as evaporative cooling relies on 
efficient elastic collisions and, even more crucially, on the ratio of the probabilities for elastic 
scattering and spin relaxation (7) which must be very large in order to prevent heating 
and trap loss. External electromagnetic fields may serve to control the rate of inelastic 
collisions. Tuning close to a Feshbach resonance has proved to be an extremely fruitful 
means of controlling atom-atom collisions 2l|] . Interestingly, it has been recently shown 22] 
that inelastic collision rates in atom-molecule collisions can be tremendously reduced in the 
vicinity of a Feshbach resonance controlled by an electric or magnetic field. 

While a large amount of work has been carried out for atom-atom and atom-molecule 
collisions, studies of molecule-molecule collisions in external fields are still scarce. Most 
clues about these more complex systems have come from atom-molecule studies. Krems and 
Dalgarno {23I1 found that the main mechanisms of spin relaxation in collisions of molecules 
with He is given by couplings to rotationally excited states mediated by the spin-spin inter- 
action. Volpi and Bohn|2^ found that spin depolarization is suppressed when the Zeeman 
splitting between incident and final states does not exceed the height of the centrifugal bar- 
rier in the exit channel. These ideas were confirmed for ^^02('^S^) + ^^02(^S^ ) by Tscherbul 
et a/|;2^], who carried out the first accurate computational study involving two diatoms. In 
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that work, the experimentally derived potential energy surface (PES) of Aquilanti et al 26| 
was employed (Perugia PES in what follows). This coUisional system is interesting since 
oxygen has been postulated as a reliable candidate for trapping and cooling^, 28 1 and 



progress in cooling this species has been achieved recently IJ, l29 |. 



The present work builds up along these lines by the investigation of the role played by 
the PES in O2+O2 collisions in the presence of a magnetic field. It is well known that 
ultracold atom-atom collisions are very sensitive to the short range of the potential 3^ . 
However, it has been recently shown 3l|] that, in the presence of inelastic scattering (i.e. 
atom- molecule collisions), peaks in cross sections around a Feshbach resonance may become 
suppressed and hence dynamics becomes rather insensitive to the details of the potential. 
This theory is tested here for a rather anisotropic molecule-molecule system such as O2 
+ O2, using a recent ah initio PES developed by Bartolomei et a/j32|. In this potential, 
electronic correlation is included by means of a high level supermolecular method in the short 
range whereas long-range interaction coefficients have been obtained from first principles 



as well 



33j. It is worthwhile to mention that inelastic rate coefficients obtained with this 



PES have proved to be highly consistent with measurements of the evolution of rotational 
populations along supersonic expansions in the temperature range 10 < T < 34 kIs^. By 
comparing present scattering calculations with previous ones using the Perugia PES [25] and 
with some additional test modifications of the ah initio PES, the effect of the potential on 
the cold and ultracold dynamics has been assessed. 

The paper is organized as follows. In Sec. II, a summary of the theory for the scattering 
between two identical molecules is given. Details specific to the ^^02 — O2 system are 
provided in Sec. Ill and in Sec. IV, results are reported and discussed. A concluding remark 
is given in Sec. V. 



II. THEORY 



We give a summary of the theory -recently developed by Tscherbul et a/25|]- for the 
scattering of two identical rigid rotor molecules in the presence of a magnetic field. 
Diatom-diatom Jacobi coordinates are used in a space-fixed (SF) frame, including the vector 
joining the centers of mass of the molecules a and 6, R, and the intramolecular unit vectors, 
fa and ff,. Intramolecular distances are fixed to the molecular equilibrium distance, Ta = 



3 



= Tg. The Hamiltonian for the interaction can be written as 



where atomic units are used = 1), / is the orbital angular momentum, fi is the reduced 
mass and V is the interaction potential or PES. The internal Hamiltonian of the molecule 
Hc,{a = a, b) is given, within the rigid rotor approximation, byjssi 



Ha = BeUl + 2^bB ■ Sa + Tsr^Q ■ S„ + (2) 

g 

where "ha is the angular momentum associated with fa. Be is the rotational constant, /i^ is 
the Bohr magneton, B is the external magnetic field and s is the electron spin. The last two 
terms in Eqj2] correspond to the spin-rotation and spin-spin interactions, parameterized by 
•jsr and Ass, respectively. Weaker interactions such as hyperfine and magnetic dipole-dipole 



are neglected (see Ref. |28l| for discussion). 

The total wave function is expanded in a basis set of SF uncoupled and symmetry-adapted 
functions 



with 



R 

Ta>Ti,lmi 



= 7Z7—^ ^ + V^lnra)) \lmi), (4) 

(2 (1 + dr^^rj) 

\lmi) being a spherical harmonics and where Ir^) represents an uncoupled function of the a 
monomer 

\ra) = \namnj\samsj. (5) 

The basis of EqJH are a well-ordered set with > r^, that are normalized eigenfunctions of 
the operator permuting the identical molecular skeletons (P: fa ff,; fb fa] R ^ —R), 
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with eigenvalue t]. These basis functions are also eigenfunctions of spatial inversion {E*: 
fa — )■ —fa] fb — ^ —fb] R —R) with eigenvalue e = (— Since the molecules under 
study are homonuclear, Ua and Ub have the same parity so e = (— 1)'. In addition to these 
symmetries, the Hamiltonian commutes with the SF z-axis component of the total angular 
momentum, so that for a given value of this projection, M, basis functions in Eqj3] must 
satisfy 

+ + + m^^ +mi = M. (6) 

Substitution of Eq|3] into the Schrodinger equation leads to the set of close-coupled equa- 
tions for the radial coefficients: 



1 _ /(/ + !) 



+ E 



u 



Mrie 

TaTtlrrii 



(R) 



(7) 



CZJiV + Ha + Hb}\<P^,^)u^,^,{R), 



where E is the total energy. It must be pointed out that the asymptotic Hamiltonian Ha + Hb 
is not diagonal in the basis (b^l^n, , due to the spin-rotation and spin-spin terms, and matrix 



elements of these terms are given in Eqs.l4 and 16 of Ref. 23|, respectively. On the other 
hand, potential matrix elements are given as a sum of a direct and an exchange coupling 



terms 



TaTbl' 



' |T/| j^Af»;e \ 



X [{Tanlmi\V\T'aTll'm[) + 7]e{Tanlmi\V\Tl^Tal'm'i)] . 



(8) 



The interaction potential depends on the total spin resulting from the coupling of the Sa 
Sb = I spins of the molecules, S* = 0, 1, 2, and can be represented as[36|: 



2 s 

V{R,fa,fb) = J2 Vs{R,fa,fb)\SM,){SM, 

S=0 Ms=-S 



(9) 



where Ms is the projection of the total spin, Ms = rris^ + m^^^. We use this representation 



in order to include directly the singlet, triplet and quintet ab initio PESs of Ref. 



32| (an 



alternative approach was followed in Ref. 25|] since the Perugia PES is given as a sum of a 
spin-independent and a spin-dependent contribution j26|). In this way, matrix elements of 
Eqj8]can be further developed as 



{TaUlmilVKrll'm'i) = 5ms,m's I](25 + 1) 



5=0 
1 1 



X 



11 S \ I I I S 

m-sa "1^6 -Ms 
X {riamn, nhUin, lmi\Vs\ n'^rUn'^ n'f^rrin' I' mi'), 



-M's 



(10) 



where (:::) are 3-j symbo 
is given in Eq.l8 of Ref. 



s. An explicit expression for {namn^n}jmnJ.Tni\Vs\n'^mniji'^pnn' I' rnii) 



25| 



Close-coupled equations (Eq|7]) are solved by means a log-derivative method j37|, |38| and 
using the basis set of EqJH in which, as mentioned above, the asymptotic Hamiltonian is 
not diagonal. At the point of imposing scattering boundary conditions and thus, obtaining 



the scattering S-matrix, it is necessary to transform to a new basis set ip? c im giving the 



eigenstates of the fragments. For each /, m; block: 



Ha + Hb 



I Mrje / I \ I ivirje 



Mri<E 



(11) 



where is the Zeeman fine structure energy level of molecule a. An unitary transformation 
of the log-derivative matrix onto the new basis is performed at the end of the propagation, 
and then scattering S-matrices and transition T-matrices are obtained in a standard way 25j . 
The integral cross section for a transition (aCb ~^ CCfe within a given (M, rj, e) block is finally 
given as 



a 



Mrje 



^(1 + ^Ca.cJ ^ ^j^Mve 



(12) 



where T is the transition matrix and /c^^^^/(2yu) = E — ec^^ +^C6 translational energy of 

the initial channel. In obtaining EqJT2| integration of the differential cross section has been 
restricted over half-space for final states satisfying = Q (see Ref. |25|). This is equivalent 
to dividing by two the cross sections integrated over full-space to avoid double counting 



when the state of the outgoing molecules is the same 
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40| 



III. COMPUTATION DETAILS 



The asymptotic Hamiltonian of Eql2] is parameterized for ^^02 by means of accurate 
spectroscopic constants 41^ : -Be=l-353 cm~^, 7sr=-0. 00396 cm~^ and Ass=l-985 cm~^. The 
three lowest states of the n= even manifold are given in Table [T] for a typical value of the 
magnetic field. Dependence with magnetic field of the combined \CaXb) asymptotic states 
is depicted in FigHJ In this work, we focus on the initial state \CaXb) = |3,3), i.e., both 
molecules are, prior to interaction, in their lowest Ifs state. Elastic and inelastic integral 
cross sections are obtained for translational energies ranging from 10~^ to 0.05 K. As we are 
dealing with collisions between identical (composite) bosons, calculations are restricted to 



the ?7 = +1 block (the role of nuclear spin can be ignored, as explained in detail in Ref. |28|). 
Note also that to study processes involving identical internal states (EqH]), calculations are 
constrained to the e = +1 parity (only even /'s in the wavefunction expansion). 

The intermolecular interaction is given by the global ab initio PES of Bartolomei et 
a/j32l|. specifically, the one referred in that work as CC-PT2 PES. Singlet, triplet and quintet 



(S'=0,l,2) potentials are given 32| by the spherical harmonic expansion 



Vs{R,fa,h) = (47r)3/2^y^^"^^\i?)A,„,,,(i?,f„,f,), (13) 

where ^AaA^A is given as a combination of spherical harmonics and A^, Af, and A are even 
integers (due to the symmetry of the four identical nuclei). The radial coefficients ^^"'^'''^(i?) 
were obtained by means of quadratures of the supermolecular ab initio energies over the 
angular variables, obtaining a total of 29 coefficients for the quintet PES and 27 for the singlet 
and triplet ones. The PESs are extended asymptotically {R > 19 bohr) using analytical 
functions (common to the three multiplicities) based on high level ab initio calculations of 
electrostatic, dispersion and induction long range coefficients [33]]. In the following Section, 
we present a comparison with calculations using the Perugia PEsj26|, which comprises just 
four radial terms (for each multiplicity) derived from a multi-property fitting analysis. To 
give a flavor of the similarities/differences between the two PES considered, we present 
in Figl2]the dependence with the intermolecular distance of the potential matrix elements 



among the Ifs state |3, 3) and the (one spin flipping) relaxation channel |3, 
elements are relevant to the mechanisms proposed by Krems and Dalgarno 



). These matrix 



23| and by Volpi 
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and Bohn2J]. Note that for initial states approaching in an s wave, conservation of M 
forbids s waves in the spin relaxation channels (see EqJH] and Ref. [24]). It can be seen 
that there are some quantitative differences in the coupling as well as in the long range 
behavior. A comparison of properties related to the van der Waals (vdW) coefficient C^^^ 
is summarized in Table [Tll 

Cross sections are computed using the code developed by Tscherbul et a/[25|, modifie d by 
us to include the hybrid log-derivative/Airy propagator of Alexander and Manolopoulos 



Related routines were taken from the MOLSCAT code 421. In this way, the log-derivative 
propagator of Manolopoulos |37| is used in the strongly coupled region (from 4.5 to 40.8 oq) 
with a fixed short step (0.04 ao), whereas the Airy propagator of Ref.[38| is used for the long 
range region (from 40.8 oq to 202. oq) with a variable step size (the ratio between adjacent 
step sizes being of 1.05). Comparing with the original code of Tscherbul et al, where only 
the log-derivative propagator was used, we found that the errors are less than 0.5% while 
the new propagation is about 10 times faster due to the smaller number of integration steps 
as well as the use of the computationally less expensive Airy propagator. 

The total wave function is expanded using a basis set comprising three rotational leve^ 



{ua, nh=0, 2, 4) and four partial waves (/=0, 2, 4, 6), equal to that employed in Ref. |25|. 



Although exact positions of the resonances might change with an increase of the basis size, 
this basis is sufficient to retrieve the main features of the collision dynamics. Regarding the 
convergence of the cross sections with the projection of the total angular momentum, M, it 
is found that for translational energies lower than 10"*^ K, just the M = 2 block calculation 
is sufficient, while for larger energies, five blocks (M =0-4) have to be summed up. For a 
single energy and magnetic field calculation, typical run times are of about 18 and 90 hours, 
respectively. 



IV. RESULTS AND DISCUSSION 

We present first the results concerning the i?-field dependence at very low energies and, 
in a subsequent section, we report those related to the translational energy dependence, 
including the transition from the ultracold to the cold regimes. 
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A. Magnetic field dependence at 1 /i K 



The magnetic-field dependence of tlie ab initio and Perugia cross sections for the Ifs state 
|3,3) at 1 is summarized in Figj3] (panels a and b). In FigJSlc we report the elastic- 
to-inelastic ratio, 7, more specifically, the ratio between the elastic cross section and those 
inelastic ones leading to untrapped states: \Q, Q) = |3, 1), |2, 2), |2, 1) and |1, 1). Note that 
new calculations with the Perugia PES were performed using the same basis set as with the 
ab initio PES (there are some quantitative changes between present calculations and those 
given in Fig.3 of Ref. 25|] were a smaller basis was employed). There are various noticeable 
differences between the two PESs. On the one hand, ab initio elastic and inelastic cross 
sections (Figs|31a,b) are much larger than the Perugia ones and, in addition, they exhibit 
more marked Feshbach resonance structures. On the other hand, although there are large 
variations of the elastic-to-inelastic ratio with the magnetic field, it can be seen that both 
PESs produce values which, on average, are of the same order of magnitude. The cases of 
very low fields {B < 50 G), where 7 is much larger for the Perugia PES, and around 1000 
G, where the ab initio value becomes very large, are discussed in more detail below. 

We discuss first the background behavior of the cross sections of FigJSl The elastic 
cross sections correspond to a background scattering length, atg, of about 118 and 32 oq 
(in absolute value), for the ab initio and Perugia PESs, respectively. These quantities are 
larger than the scattering lenghts purely due to the vdW potential js^, a, of 22 and 24 
ao, respectively. The particularly large value of the ab initio elastic cross section can be 
explained by existence of a close quasibound state varying with magnetic field at the same 
rate than the entrance channel. Regarding inelastic cross sections, the ab initio one is on 
average about 10 times larger than the Perugia result. This difference can be qualitatively 



rationalized by resorting to the analytic van der Waals theory 



43 



44j |. which takes the 



solutions of the vdW potential 45| as the reference for the multichannel quantum defect 



theory 46| . A key parameter in that approach is the short range squared amplitude of the 
entrance channel wave function, which near threshold is proportional to|43| 



limCbgiko) 



1+(1-^ 

a 



(14) 



ko being the wavenumber of the incoming channel. Since inelastic cross sections are pro- 
portional to C(,g(A;o)~^|46|, EqJT^ implies that the value of abg affects the threshold behavior 
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of the inelastic cross sections. It follows, then, that the very large ab initio inelastic cross 
sections are explained by the magnitude of the corresponding background scattering length. 
Within this framework, one can expect that the elastic-to-inelastic ratio becomes less sensi- 
tive to abg than the cross sections themselves, since both elastic and inelastic cross sections 
are approximately proportional to a^g. This is the result of FigJSlc, where the average value 
of 7 is about the same for both potentials. 

We now turn to discuss the resonant structures of Fig 131 At this point, it is convenient to 
mention the work of Hutson 3l| who analyzed the threshold behavior of Feshbach resonances 
in the presence of inelastic scattering. He found that -in contrast to the case of a pure elastic 
scattering- resonance peaks may be significantly suppressed and, in this way, the coUisional 
process may become insensitive to the details of the potential. With this in mind, the 
profiles obtained in Figl31 are rather unexpected given the considerable anisotropy of the 
O2 — O2 interaction. In connection with this issue, let us digress for a while and study the 
resonance patterns for a purely elastic scattering event, as is the case of the magnetic field 
dependence of the lowest high field seeking (hfs) state |1, 1) (see FigUl). The result for the ab 



initio PES at 1 /iK, using a a reduced 
directly compared with Fig.4 of Ref. 



_)asis {nmax=4:,lniax=4:) , is shown in FiglHand can be 
25| . For both PESs, a high density of very pronounced 
resonances is obtained. For the ab initio PES there is a slightly larger number of peaks, 
and some of them are wider. Also, the baseline of the ab initio cross section is much larger 
than the Perugia one, as occurs for the Ifs state. A similar density of quasibound states 
is expected when the entrance channel is the Ifs state, but presence of inelastic channels 



substantially modify the resonance lineshapes 



31^ 



To show this, it is convenient to write 



down the behavior of the S matrix in the neighborhood of an isolated resonance |31l. |4 

where k and j are the incoming and outgoing channels, respectively, S*^^ is the background 
S matrix, E is the total energy, is the resonance position, F^; is the resonance width and 
(complex) gEi involve couplings between resonance and eM,.e. . wavefunctionfl, such 
that the partial width for channel i is given d&TEi = \gEi\^ and F^; = F^;^. A key point in 

Hutson's argument is that gEk elements are proportional to the square root of the incoming 

1/2 

channel wavenumber k^ . Then, as ko decreases and if the resonant state is also coupled 
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to inelastic channels, the radius of the circle described by 5',^ drops to zero and peaks in 



cross sections become significantly supressed 



3ll ]. The analytical vdW theory gives a more 



detailed threshold behavior of the qeu elements, as they become proportional to the square 
root of EqJT^ Hence, if abg is sufficiently large, qeu will tend its threshold value (zero) 
rather slowly, and as a consequence, more pronounced peaks in the cross sections can be 
obtained. This explains why we find a marked resonance structure, especially for the ah 



initio PES. Nevertheless, as noted in Ref. 



3l| . a relatively large ratio between elastic and 



inelastic partial widths is also needed in order to obtain pronounced resonance profiles. It 
is reasonable to expect that, among all the quasibound states that should be crossing the 
Ifs state, only some of them will have particularly large elastic partial widths, so only a few 
marked resonance features will "survive", as in fact it occurs (FigjS]). 

We have just seen that a large a^g enhances the short range couplings between the res- 
onance and the incoming wavefunctions. In this situation the dynamics must become very 
sensitive to the short range region of the potential. In order to study the role played by the 
short range vs. the long-range features of the intermolecular potential, we have performed a 
test calculation where the long-range anisotropy of the potential is switched off. To this end, 
the ah initio PES has been modified by imposing, for i? > 19 Oq, an exponential decay of all 
radial terms of Eq{T3]except the isotropic one (Aq, \, A) = (000). The new cross sections are 
compared with those corresponding to the correct long range behavior in FiglSl This figure 
clearly shows that the resonance structure is rather insensitive to the long-range anisotropy 
of the interaction and, therefore, short range couplings must be playing a dominant role. 

Finally, it is interesting to note from FigJSlb that, for the ah initio PES, there is a 
significant suppression of inelastic scattering for magnetic fields ranging from 750 to 1500 G. 
This feature must be related with the prominent resonance at about 600 G and it must be due 
to interferences between the background and resonant S matrices leading to asymmetric line- 
shapes of the state-to-state cross sections 49j . Note that this reduction entails a considerable 
increase of the ratio 7 for a wide range of magnetic fields. A similar behavior (with an even 
larger suppression of inelastic scattering) has been found in "^He -|- ^^02 magnetic Feshbach 



resonances 



Analogously, it is also worth mentioning that, for the ah initio results, the elastic scat- 
tering on the left- hand- side of the resonance at about B= 30 G is suppressed. This feature, 
already present in Figj3l can be more clearly seen in FiglSl where the ah initio elastic cross 
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section becomes very small around 10 G. In this case, the corresponding ratio 7 becomes 
much smaller than expected (from the well known effect of suppresion inelastic scattering 
due to centrifugal barriers |24J. l25l . l28|). 



B. Translational energy dependence 

In Figjni dependence of the cross sections with kinetic energy is given for several selected 
values of the magnetic field. In agreement with predictions based on the analytical vdW 
theory 43| , two very different regimes are noticed for energies larger or smaller than Eydw ~ 



10 mK (see Table [ITl)- For the higher energy range, elastic and inelastic cross sections 
exhibit a weak dependence with the field, the Perugia ones being larger than their ab initio 



counterparts, in consistency with previous studies at higher energies [39|]. For energies lower 
than the crossover [E^dw), cross sections become more dependent on the magnetic field. 
This is mainly due to the effect of the resonances in the ultracold regime, but in the case of 
the Perugia PES, supression of inelastic cross sections at low fields (due to the centrifugal 



barriers 



25 



28| ) also plays a role. 

It is interesting to highlight that a relatively high value of the elastic-to-inelastic ratio has 
been obtained between 1 and 10 mK in the ab initio calculation at 1000 G (Figj6lc). This 
result is related to the asymmetry of the lineshape and the suppression of spin-changing 
processes on the right-hand-side of the resonance at 600 G and 1 /i K, discussed above 

(FigED. 

A more detailed study of the ab initio cross sections for low values of the field [B < 50 
G) is given in Fig|71 An impressive dependence with B is noticed for energies just below 10 
mK. Between 1 and 10 mK, complicated resonance structures are seen which are particularly 
acute for the elastic cross section. These features are related to the prominent resonance 
around 30 G at much lower energies (reported in Figj3] and more clearly seen in Figj5]). 
In other words, they are expressions -at several different energies and magnetic fields- of 
the same quasibound state. For instance, note the resemblance between the asymmetric line 
shapes of the elastic cross section at B= 1 and 5 G and between 1 and 10 mK (FigjTla), with 
the magnetic field dependence at much lower energies for fields i? < 30 G, shown in FigjSl A 
detailed tracking of these resonances would involve non-trivial lineshape fittings and has not 
been attempted here. On the other hand, it should be noted that, for the range of magnetic 
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fields of FigJT] and up to translational energies of at least 1 mK, spin- changing collisions 



should be suppressed due to existence of centrifugal barriers in all outgoing channels 



25 



28|. 



In FigJTlb it can be seen that, except for the lowest value of B (1 G), such a suppression 



does not occur, in contrast with the results using the Perugia PES (see Ref. |25l| and FiglHlb). 
This must be due to a significant tunneling through the centrifugal barriers for energies/fields 
close to the resonance. Consequently, the ratios 7 are particularly small for this range of 
fields (FigHb). 

A further analysis of the sensitivity of the elastic-to-inelastic ratio to the details of the 
PES has been performed. We have artificially modified the anisotropy of the present ab 
initio PES by multiplying all the terms in the spherical harmonic expansion (Eq JT3|) -except 
the isotropic ones- by a factor /3 ranging form 0.98 to 1.02. In Figj8]we show the results for 
different translational energies and magnetic fields. In can be seen that, while for 20 mK 
there is not a strong variation of 7 with /3, for lower energies (1 mK and 1 fiK), this ratio 
changes tremendously with the anisotropy of the potential. In the new calculations {(3 = 
0.98 and 1.02), no nearby resonances appear for the energies/fields considered and hence, 
results are "more standard", i.e., very large values of 7 are now attained for low values of 
the field {B < 50 G), in agreement with the expected suppression of inelastic scattering, 
but smaller 7's are obtained for B= 1000 G. However, note that, contrarily to a first order 
perturbation theory, the largest ratios are obtained with the most anisotropic PES (/3= 
1.02). 

V. CONCLUDING DISCUSSION 

We have performed a detailed study of cold and ultracold molecule-molecule collisions in 
the presence of a magnetic field for a system with a significant anisotropy such as O2+O2. 
A thorough comparison has been made between a high quality ab initio PES and previous 



studies 251] where a a different PES was used. Several interesting findings have emerged 
from this approach regarding the anisotropy as well as the relative influence of long and 
short components of the interaction. For the ab initio PES, a large background scattering 
length gives rise to pronounced resonance structures in the ultracold regime (translational 
energies < 10 mK). As a consequence, the ratio between elastic and inelastic cross sections, 
7, is very dependent on the magnetic field as well as on the short range anisotropy of the 
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PES. Therefore, quantitative predictions for this important parameter become rather risky. 
However and as a general trend, we can indicate that high values of 7 could be achieved 
in the vicinity of asymmetric Fano resonances, or for low fields, -B < 50 G. Note that the 
maximum temperature that can be held in a trap with such a depth would be of about 1 



mK 



25 



5Q|. 



A key issue is the large density of quasibound states of the O2+O2 system, best illustrated 
in the magnetic field dependence of the elastic cross sections of the lowest high field seeking 
state. In view of this, having obtained a large background scattering length does not seem a 
rare event. Present behavior might be characteristic of a range of molecule- molecule systems 
as well, that is to say, as the number of degrees of freedom increases, a larger density of 
quasibound states, including near threshold resonances, can be expected[51i]. which in turn 

53] have proposed an efficient 



makes dynamics richer. Very recently, Suleimanov and Krems 
method for locating Feshbach resonances in external fields. The new method could be very 
useful for the comparison of spectral patterns obtained from different potentials or between 
different molecular systems. 
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TABLE I: Energies (in K) and coefficients (in tlie basis of Eql5]) of tlie tliree lowest states of ^^©2 
(EqM) for a magnetic field B= 100 G. 
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TABLE II: Parameters associated to the long range behavior of the ah initio and Perugia potentials: 
isotropic vdW coefficient Cg^", scale length {Rvdw) and energy {E^dw) of the analytical vdW 



theory 



431 ] ■ and height of the d-wave centrifugal barrier. Bmin is the critical magnetic field for 



which the |3, 3) — |3, 1) Zeeman splitting becomes larger than the d-wave barrier. 

ab initio Perugia 
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FIG. 1: Internal energies of ^^02(Ca) +^''' 02{Cb) as functions of magnetic field. In this work, 
molecules are considered to be initially in their Ifs states \Ca, Cb) = |3, 3). Open and closed symbols 
indicate critical values of the field for which d and g barriers, respectively, become open for the 
different outgoing channels. Note also that calculations of FigS] refer to the hfs state |1, 1). 




FIG. 2: Comparison of the ab initio and Perugia potential matrix elements among the fragment 
states |3,3) and |3, 1) for a magnetic field 5=100 G. The long range behavior is compared in the 
inset. Note that orbital angular momentum for the entrance (|3,3)) and outgoing(|3, l))channels 
are and 2, respectively. 
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FIG. 3: Magnetic-field dependence for collisions of ^^02 +^'' O2 in the initial Ifs state |3,3) and 
translational energy of 1 /xK. a) Elastic cross sections; b) Total inelastic cross sections and c) Ratio 
7 between elastic and inelastic (untrapping) cross sections. Blue and red colors correspond to using 



the ab initio 



34] and the Perugia 



261] PESs, respectively. 
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FIG. 4: Total cross section vs. magnetic-field for the hfs state |1,1) at a translational energy of 
1 /uK and using the ab initio PES. Note that only the elastic channel is open. The result can be 



compared with Fig. 4 of Ref. 



25l |. corresponding to the Perugia PES. 
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FIG. 5: Effect of the long-range anisotropy of the ab initio PES: Magnetic-field dependence of 
elastic and total inelastic cross sections for the Ifs state |3, 3) at lOuK. Thick lines joined by filled 



33l |. while dashed lines joined by 



squares show results using the correct long range anisotropy 
open squares correspond to calculations where the long range anisotropy of the interaction has 
been switched off. 
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FIG. 6: Translational energy dependence of the collisional processes for the Ifs state |3,3): a 
comparison between ab initio and Perugia PESs for different values of magnetic field, a) Elastic 
cross sections; b) Total inelastic cross sections and c) Elastic-to-inelastic ratio 7. 
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FIG. 8: Study of the sensitivity of the elastic-to-inelastic ratio 7 to the anisotropy of the ab initio 
PES. All terms of the spherical harmonic expansion are multiplied by /3 except the isotropic one. 
The effect is shown for different translational energies and magnetic fields. 
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